Universal interface width distributions at the depinning threshold 
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We compute the probability distribution of the interface width at the depinning threshold, using recent power- 
ful algorithms. It confirms the universality classes found previously. In all cases, the distribution is surprisingly 
well approximated by a generalized Gaussian theory of independant modes which decay with a characteristic 
propagator G(g) — ^, the roughness exponent, is computed independently. A functional renormal- 

ization analysis explains this result and allows to compute the small deviations, i.e. a universal kurtosis ratio, 
in agreement with numerics. We stress the importance of the Gaussian theory to interpret numerical data and 
experiments. 



The scaling properties of driven elastic interfaces in random 
media play an important role in a wide variety of physical 
situations, ranging from stochastic surface growth to domain 
walls in disordered magnetic materials, the spreading of fluids 
on rough substrates, and the dynamics of cracks {Hl^. These 
problems share many features with critical phenomena and 
provide a challenge for theoretical approaches to disordered 
systems and non-equilibrium phenomena 

Here we study interfaces described by a scalar height func- 
tion h{x), where x is the d-dimensional internal coordinate. 
We measure the deviation from the mean position as u(x) = 
/i(x) — (h), where (. . .) stands for the spatial average over all 
X of a given interface (cf. Fig.[0. The mean square width of 



<h> 




a single interface, w ({u(x)}) 



can be used to char- 



acterize its roughness, and explore universal properties: After 
averaging over the ensemble of interfaces, w'^ grows with the 
lateral extension L of the system as 

~ cx i^C for ^ 



oo 



(1) 



where C is the roughness exponent. 

An interesting property is that, for positive (, uP' fluctu- 
ates even in the thermodynamic limit |9, 10, 11]. This means 
that the long-range geometric features of the interface are not 
characterized by the roughness exponent alone, but require 
the complete probability distribution P{w'^). P{w'^) has been 
computed for several linear stochastic growth equations with- 
out disorder as the Edwards-Wilkinson model, the Mullins- 
Herring model, and the 1-d KPZ model lliol ll ill . In these 
models, the probability distribution P{w^) can be rescaled 
into a form independent of system size and of microscopic 
details 



P(u)2) = {l/w^)<i>{w^/w^) for L 



(2) 



Although w"^ may contain a non-universal scale, the function 
^{z) is universal. It has been argued that the shape of P{w'^) 
can thus be used as a sensitive tool, distinct from to distin- 
guish between different universality classes |9, IJiJ 
Furthermore, $(2;) is expected to converge to a 5-function 
above the upper critical dimension due- This has motivated at- 
tempts to determine due for e.g. the KPZ equation 1 13]. Proba- 
bility distributions of order parameters have received much at- 
tention for related models such as polymers, spin glasses, and 



FIG. 1: Example of a 1 + 1-dimensional periodic interface h{x) 
(random walk) with mean value (ft), and u{x) — h{x) — {h). 



random diffusion ] 14]. The quantity we study here, P{w'^), is 
the distribution of the lowest order observable which tests the 
whole function h{x) for < a; < L. It appears as a funda- 
mental quantity in disordered systems. 

The aim of this Letter is to compute the width distribution 
(WD) $(z) for elastic interfaces driven in random media, ex- 
actly at the depinning threshold, numerically and from field 
theory. As in the linear problems treated earlier, we confirm 
the existence of universal properties in various dimensions d 
and with several functional forms of the elasticity. The sur- 
prising finding is that in all cases P{w'^) (i.e. its shape) is ex- 
tremely well approximated by a simple generalized Gaussian 
approximation (GA), without any fit parameter, and depends 
only on (, which is determined independently. This suggests 
that the complicated morphology of interfaces (cracks, do- 
main walls, etc.) may be rendered by a simple ansatz of inde- 
pendent modes with a characteristic decay. This may have im- 
portant consequences for the analysis of numerical and experi- 
mental data. Our numerical results are then understood within 
a functional renormalization group calculation, detailed in a 
companion paper L15il . 

We consider the zero temperature equation of motion of an 
interface given by 
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(3) 



where the functional E({h, x}) represents the total energy in- 
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corporating potential energy due to the driving force /, the 
short-range correlated disorder force 77(x, h), as well as its in- 
ternal convex elastic energy Eg]. Equation is non-linear, 
and has not been solved exactly. We are interested in the de- 
pinning limit (/ = /c) where the velocity of the elastic man- 
ifold goes to zero. We use periodic boundary conditions, and 
recall that the WD $(2), although independent of small scale 
details, does depend on the boundary condition at large scale. 
For our numerical study we use a very efficient algorithm 
which directly determines critical forces as 
well as the critical interface /ic(x) for a wide range of models. 
In particular we calculate the WD for interfaces of dimensions 
d — 1 and d = 2, where the elastic energy has the harmonic 
form Ec\{{h,x.}) ^ (V/i)^. We have also tested the univer- 
sality of ^{z) in d = 1, by means of a directed polymer model 
with an anharmonic quartic elasticity, and for a lattice model 
with hard local constraint, which have the same C = 0.63 fl^l- 
As expected, ^{z) is always size independent and the WD as- 
sociated to non-harmonic models can be distinguished from 
the one resulting from an harmonic elasticity. The harmonic 
models, in fact, have an exponent ^ — 1.2, and thus belong to 
a different universality class. 

For our field theory calculation we use the functional renor- 
malization group method (FRG) originally developped to one 
loop to describe the model with harmonic elasticity and cor- 
rect the predictions of dimensional reduction R*^. Recently 
a renormalized field theory was constructed to 2-loop order 
0] which overcomes the deficiencies of the 1-loop analysis; 
notably to distinguish between statics and driven dynamics, 
and to account for the large values of the roughness exponent 
^ measured e.g. in as compared to an earlier con- 

jecture 1 5] C = (4 - d)/3. We find that the FRG both suggests 
the OA as a lowest order approximation in e = 4 — d and al- 
lows to define and compute universal ratios which probe high 
cumulants of P(w^) and deviations from the GA, and are thus 
more sensitive to details of the universality class. The simplest 
of them is the generaUzed kurtosis 
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where the subscript c indicates the connected expectation 
value. R is found to be small but non-zero. This directly 
proves that the correct description of interfaces must go be- 
yond the independent-mode picture. 

To introduce the Gaussian approximation in the most ele- 
mentary way, we first recall [9] the simple periodic random 
walk of size L, with a Fourier expansion 

u[x) = 2_^anCOS [^—x j + bnSm [—j—x j . (5) 

The standard Gaussian probability measure associated with u. 



i.e. Vlu] oc exp[-i dx{^^)'^] gives 

(7rn)2 



V[u\ oc exp 



E 



The probability distribution P{w'^) 

is obtained from the generating function of its moments 
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Writing Eq. (|8ji as an integral over a„ and 6„, we obtain 



Z{X) = n / da„d5„e"^(«"+''")e-t(-?.+''^) 



(6) 



(7) 



(8) 



W{X) = 
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Z{0) 
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(9) 



For the random walk Eq. (|6}, P{w^) can be obtained exactly 
by inverse Laplace transform of Eq. (|9}: 



P(w^) 



(10) 



n>0 



Using vP' — — ^^|a=o = "^l- * I10> can be written in a 
scaling form 

$(z) = 'uP-P{w^) , z = v? j-uP- 
3 

n>0 

The size dependence thus appears only through the average 
width . We can generalize Eq. (|6j, where each mode a„, 6„ 
has a weight cx n-^, to an arbitrary function of independent 
Fourier modes 



T'gausslw] OC exp 



n>0 ^ ' 



(12) 



which, in real space, corresponds to 



"Pgaussi"] (X exp 
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AxAyu{x)G u{y) 



(13) 



The function Gxy — u{x)u{y) is the exact disorder-averaged 
2-point function and can be computed from numerical data. 
This allows to obtain T'gauss even for a finite system. In 
the thermodynamic limit, T'gaussiw] is obtained from the be- 
havior of Gxy = Gx-y for large \x — y\ (small q), where 
G{q) ~ Cjq^^'^'^. This means that a single observable, 
fixes T'gaussiw] on large scales. 
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FIG. 2: Scaling function $(2:) and $gauss(2) for: (1 + l)-d harmonic 
(L = 256, 2 X 10^ samples, C = 1.25)~left; (1 + l)-d an-harmonic 
(L = 256, 2 X 10^ samples, ( = 0.63)— middle; (2 + l)-d harmonic 
(L = 32, IC samples, — 0.75) — right. The scatter in the numeri- 
cal data is mostly due to binning. Notice that, for d = 1, the typical 
value of z is much smaller than its average 1=1. 



We again determine the generating function for the mo- 
ments, but for arbitrary ( and d: 



W{X) 



n 



1 + 2AG(g) 



-1/2 



(14) 



where A = \/L, q = 2Trn/L, n E l/. Due to the symmetry 
g <^ — (?, no fractional power appears in Eq. ( I14> . as in Eq. (|9jl, 
where the exponent — 1 stems from the double sum over the a„ 
and 6„. An explicit sum over poles allows to obtain $gauss(z) 
for all C and d with excellent precision. All GA interfaces 
\u(x)] can be directly sampled by Monte Carlo methods. For 
details, including the extension to open boundary conditions, 
see 1 12. 20] 

In Fig. 121 we compare, for different models, the exact scal- 
ing function (f>(z) to $gauss(z), using G{n) = Cjvf'^'^^. The 
roughness exponent was previously obtained using both field 
theory |7] and numerical methods |17|. The agreement be- 
tween $ and $gau,ss is clearly spectacular The scatter of the 
data, visible in Fig.|2l is mostly due to the finite width of his- 
togram bins. 

Tiny — yet significant — differences between $ and <i>gaus,s 
are best resolved in the integrated probability distributions, 
which need no discretization. The difference between the in- 
tegrated distributions is 
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(15) 



where H{z) = dt^(t) is the fraction of samples with a 
renormalized width below z. In Fig.|3l we show AH{z) ob- 
tained from N = 2 10^ independent samples. Statistical fluc- 
tuations in this quantity are of order and the signal- 
to-noise ratio would be smaller than 1 if iV was an order of 



FIG. 3: Difference between the integrated distribution functions of 
$ and <I>gauss (with ( = 1.25) obtained from 2 x 10^ independent 
interfaces at L = 256 (continuum line) and L 
the d = 1 harmonic model. 



64 (dashed line), in 



magnitude smaller The absence of systematic finite-size ef- 
fects shows that the asymptotic regime of large interfaces has 
been reached and thus to conclude that the exact distribution 
for large systems is not Gaussian. 

We now discuss the field theoretical calculation. To low- 
est order in perturbation theory, we show that the generalized 
Gaussian approximation appears naturally. This is instructive 
since it identifies the diagrams which are obtained by assum- 
ing the theory to be Gaussian, albeit non-trivial, since it in- 
volves a non-trivial roughness exponent C. Using dynamical 
field theoretic methods f?"], one starts again from the Laplace 
transform T4^(A) and expands in powers of the correlator of 
the pinning force A(u). To lowest order one finds il^ that 
logW{X) is the sum of all connected 1-loop diagrams. The 
loop with N disorder vertices and insertions of is 



2N 



E 



-2AA(0) \ 



N 



(16) 



where the sums over q thus run over a d-dimensional lat- 
tice with spacing and the 0-mode is excluded, as ap- 
propriate for periodic BC. Resumming ( II6I 1 over q would 
give Eq. il4\ with G{q) ^ l/'z"', i-e. the dimensional re- 
duction (Larkin) result. In fact, the ERG tells us that the 
calculation should be performed with the running disorder 
A(0) Ai(0) = e(^-2C)'A*(0) where A*(0) is the (non- 
universal) value of the fixed point |7]. Eor the present case of 
periodic boundary conditions and momentum infrared cutoff, 
one can replace I log(l/g), and finally obtains Eq. il4l 
with G{q) = C/g'*+^^. This calculation is valid to domi- 
nant order in e = 4 — d, i.e. near d = 4. If the same class 
of diagrams are resummed in any d it leads to the GA, as we 
now illustrate considering e.g. the second connected cumulant 
of the WD. This cumulant is not connected w.r.t. h, and thus 
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there is an exact relation: 

(^" = (^-(^^ = 2(1 + 7?)/ Gly. (17) 

J x,y 

The first term results from Wick's theorem and would be the 
full result if the measure were Gaussian. Analogous formulae 
exist for higher cumulants, and if the measure of h is purely 
Gaussian can be resummed into Eq. Even though the 

GA is not exact, the deviations, given by the last term in ( I17> 
are expected to be small; indeed they are of order e*. Thus the 
GA is already exact to the two lowest leading orders and e"^, 
which explains why it is so accurate even in low dimension. 

The calculation of the deviations using the field theory is 
delicate The kurtosis R in Eq. @ which character- 

izes the importance of non-Gaussian effects is found to be 
R = — 0.13e^ to lowest order for small e — 4 — d. It is easy 
to see that this strongly overestimates R in low dimensions. 
Another method is to work in fixed dimension and to truncate 
to one loop, yielding R = -0.036 (d = 3), i? = -0.048 
(d — 2), R ^ —0.01 (d — 1), which in view of the numeri- 
cal results below seems to underestimate R. The small values 
obtained in low dimensions arise from kinematic constraints 
in the diagrams, presumably a genuine effect indicating large 
corrections from higher orders in e to the above O(e^) even 
in d = 3. Note that the sign of the result indicates a distri- 
bution more peaked than a Gaussian and is in agreement with 
the only other known exact result [21,1 for the (random field) 
statics in d = 0, i? = -0.080865... 

We have computed from Eq. (0J the generalized kurtosis 
function, in a model-independent way. We have checked on 
the 1 — d harmonic model that, as ^{z), i? is not affected by 
finite size effects and - using 10^ samples - we find: R = 
-0.054 ± 0.002 (1 - d harmonic L = 256); R = -0.067 ± 
0.002 (1 - d anharmonic L = 64); R = -0.053±0.002 (2-d 
harmonic L ~ 32). As proven by ERG calculations, R is 
small but definitely different from zero. Direct information on 
the non-Gaussian effects can also be obtained from the Fourier 
transforms of the interfaces u{x) in Eq. (|5}- In Fourier space, 
for d = 1, the expression of R is: 

We remark that R detects correlations in the disorder-averaged 
fourth moments |u((7i)p|u(g2)p, which cannot be expressed 
simply through the second moments |u((7)p. 

To summarize, we have computed both numerically and 
within field theory the width distribution of critical configu- 
rations at depinning, with consistent results. The shapes of 



the distributions are strongly dominated by the value of On 
the other hand, it will be difficult to distinguish different uni- 
versality classes from the forms of ^{z), if their roughness 
exponents are similar Other universal quantities, such as R 
defined here, directly involve the non-Gaussian part of the dis- 
tribution. Their precise determination however requires more 
work, both numerically and within field theory. Also, since 
the WD is so tied up to C,, finite size effects in both quantities 
are connected. Finite-size effects will need to be well under- 
stood in order to resolve open issues 0, [l^ E2I1 concerning 
due for the anisotropic depinning class. It would be interest- 
ing to carry similar calculations on other pure, and disordered 
models, in particular for the equivalent static system. 
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